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ABSTRACT 


A boundary point collocation procedure is applied 
to find stress intensity factors of single edge crack 
problems and centre crack problems with traction and/or 
displacement prescribed boundary conditions. For solv- 
ing the crack problems we made use of William/ s stress 
function [4] which gives the solution in terms of tri- 
go name trie terms with undetermined coefficients. The 
boundary conditions are derived from William's stress 
function [4] by successive partial differentiation. 

Thus in our present work we make use of direct boundary 
conditions, for collocation procedure. The usage of 
William's stress function prescribed for edge cracks is 
extended to centre crack problems as well in our work. 

The results thus, obtained are compared with the exist- 
ing. results and good agreement is observed. A butt weld 
improper penetration problem is idealized as a centre 
crack problem and its stress intensity factor is also 
calculated. 



CHAPTER - 1 


INTRODUCTION 

1.1 General Introductions 

Fracture mechanics is a discipline of engineer- 
ing dealing with the analysis and prediction of mecha- 
nical and structural behaviour of a cracked structure 
when it is loaded. Many structures may be made essen- 
tially csrack (or flaw) free. However/ the service con- 
ditions for many complex and high performance structures 
(air planes, offshore drilling rigs, bridges etc.) in- 
volve repeated loading, which often limits the life of 
a structure. The finite life is a result of repeated 
loads (and/or temperature fluctuations) which induce 
tiny cracks formed by the process called fatigue, often 
at some geometric discontinuity (for example, a hole). 

In service those tiny cracks grow leading to a cracked 
structure. 

Many structures are typically fabricated by 
welding. Common weld designs often introduce geometries 
that are effective cracks. For example, in butt-welding 
of thick plates due to improper penetration, we g,et 
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a structure with an effective crack at the centre. 

Further the process of welding in itself, inevitably 
leads to a compromise between slow welding speed which 
minimizes the disturbance of the metal, and a high weld- 
ing speed to increase productivity. The result is a 
structure that is unintentionally cracked by the time 
that it enters service. 

Past experience shows that whenever unexpected 
failures (at design conditions) have occurred they always 
involve a flaw or a crack. Since it may not be possible 
to prevent cracks occurring in service it is necessary 
to have a means of assessing their effect under working 
condition in a given structure. It can well be stated 
here that a cracked structure " need not necessarily' ' 
have failed. Certainly it is true that some cracked 
structures have failed for all practical purposes. However, 
it is equally true that many large complex structures 
contain 'cracks' and inspite of them will provide satis* 
factory, safe and economical performance for many years 
to come. There is no implication here that cracks or 
flaws are desirable or intentionally introduced into 
structures. There is only the recognition that many 
structures are unintentionally flawed and we must deal 
with this reality. This justifies the investment of 
money, time and effort on the study of cracked structures. 
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Fracture mechanics provides a basis for quanti- 
fying the behaviour of cracks during both the propoga- 
tion phase and the final failure. The important para- 
meter in this field of science is the " Stress Intensity 
Factor"/ / K / , which is a measure of the magnitude of the 
stress occurring in the immediate vicinity of the tip of 
a crack. The stress intensity factor is a function of 
the loading on the cracked configuration and of the size 
and shape of the crack and other geometrical boundaries. 
Thus, there is a need to know, K, for cracked configura- 
tions which might be encountered in practice. 

1.2 Literature Survey; 

Exact mathematical solutions for stress inten- 
sity factors have been determined by sophisticated mathe- 
matical techniques, like perturbation techniques and the 
method of integral transforms, for a number of idealised 
geometric configurations and loading conditions. But 
in most complex engineering situations where it is diffi- 
cult to obtain theoretical solutions. Various numerical 
methods must be used to obtain approximate solution. 

The two most general numerical methods for the elastic 
analysis of cracked solids are boundary collocation pro- 
cedure and the finite element method. Finite element 
method is the most general of the two in that it can be 
applied to curved cracks and comples geometries, and 
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loading conditions. Collocation procedure can be applied 
to straight cracks and fairly simple external boundaries. 
However, for classes where a boundary collocation pro- 
cedure is applicable, it offers the distinct advantage 
over the finite element method of greater accuracy, less 
input preparation and a relatively smaller system of 
simultaneous equationsto be solved. 

The boundary collocation procedure can be for- 
mally defined as follows : using a stress function, 
governing the stress field of a cracked solid, the boun- 
dary conditions are satisfied at a finite number of 
discrete boundary points on the outer boundary, of the 
cracked body. This is done by truncating the terms of 
the stress function such that the number of arbitrary 
constants of the stress function is equal to the number 
of discrete boundary conditions. The values of the 
constants are chosen such that the boundary conditions 
are satisfied exactly at the discrete boundary points. 

W.K. Wilson [l] solved a plane elasticity problem 
of a finite rectangular plate having a centre crack with 
arbitrary orientation under tensile loading, by means 
of an Airy's stress function. He reduced the boundary 
conditions into one of the three following forms: 

(a) surface tractions 
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(b) the normal and tangential derivative of the 
stress function of the boundary points, 

(c) the Airy's stress function itself, and its 
normal derivative of the boundary points. 

He observed that boundary collocation method using the 
boundary conditions in (c) gives the most accurate 
results. He also claims that stability, simplicity and 
relative independence of results w.r.t. the position 
of boundary points are the method's main features. 

M. Isida [2] studied a centre crack problem with 
the crack exactly perpendicular to the loading conditions 
using a complex stress function. J.E. Snawley [3] 
analysed an edge crack problem in tensile loading using 
a real stress function formulated by M.L. Williams [4]. 
Both Isida and Srawley used the boundary conditions 
for collocation procedure in the form as given in (c) 
and they obtained very good results. 

1.3 Objective and Scope of Present Works 

It may well be stated here that when we have 
both the displacement and traction prescribed problems 
over the boundary it will be convenient to use them 
directly as they are> with some special technique to 
achieve stability and accuracy even at the, expense of 
relatively higher computational time. This is not 
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unjustifiable considering the advent of computational 
machines in the present days. 

An attempt is made in our present work to use 
the boundary conditions involving second order partial 
differentiation of the William's stress function, for 
both edge crack and centre crack problems. By this we 
incorporate the traction and displacement boundary con- 
ditions directly in our problems. The centre crack 
problem is analysed treating it as two edge crack pro- 
blems. Also a butt weld improper penetration problem 
is idealized as a centre crack problem and is analysed 
using the boundary collocation procedure. 

In the chapter [2] we describe the formulation 
of the problems with the necessary boundary conditions. 
The required equations in this connection are derived. 
In chapter [3] we present the method of solution along 
with results and discussions. 
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CHAPTER - 2 

CRACKS UNDER TENSILE LOADING 

2.1 Introductions 

In this chapter William's stress function [ 4 ] 
governing the stress field around the crack tip of an 

1 

edge crack is reviewed, along with the general stress 
field around a crack tip and the modes of deformation. 
The relationship betx^een the stress intensity factor and 
the first coefficient of William's stress function is 
derived. William's stress function is used to analyse 
edge crack problems with traction and displacement boun~ 
dary conditions, A centre crack problem Is analysed 
treating it as two edge crack problems. A butt weld 
problem with improper penetration is treated a centre 
crack problem and the necessary equation in this regard 
are also derived. 

2.2 Stress Field Near a Crack Tip; 

The linear elastic solution for the stress 

field around a crack shows that the stress components 

-7. . at a distance 'r' from the crack tip can be 
1 3 1 s 

expressed in the following form (fig. 2.1.a) 
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* -j -« = fj j (©) + 'other terms' (2.1) 

J y 2 TT r 

where K, is the stress intensity factor and (r, 9) are 
the polar co-ordinates. If the point (r, 9<) 'is suffi- 
ciently close to the tip, i.e., r << crack length, the 
'other terms' in the above equation are negligible com- 
pared to the first term. The above elastic solution 
predicts infinite stresses at the crack tip (r = 0), 
which can not occur in practice, since plastic flow 
takes place in the highly stressed region near the 
crack tip. However, if the region of plastic flow is 
small compared to the region over which l /r~” term domi- 
nates the stress field, it may be assumed that, the 
behaviour of the crack is determined by the elastic 
stress intensity factor. This assumption forms the 
basis of linear elastic fracture mechanics. 

2.3 Modes of Deformations 

The stress intensity factor, K, is a function 
of the loading on the cracked conf iguration and the size 
and shape of the crack and other geometrical boundaries. 
It has the dimensions of stress x x / length . There are 
three distinct types of cracking modes which are charac- 
terized by different symbols for the stress intensity 
factors. These are illustrated in fig.(2.1.b). Mode I 



characterized by is Known as the opening mode. Mode 
II characterized by is known as the sliding mode and 

mode III characterized by ^ XI1 is known as the tearing 
mode. The stress intensity factor for each mode can be 
formally defined as follows s 
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( 2 . 2 ) 

out of these three, the opening mode, is the most common 
mode found in practical situations and in our present 
work we restrict ourselves to determine the stress inten- 
sity for the opening mode. 

2.4 William's Stress Function; 

In his paper M.L. Williams [43 has formulated 
a stress function in polar co-ordinates governing the 
stress field around the tip of an edge crack. Besides 
satisfying the biharmonic equation the function also 
satisfies the stress free conditions along the crack 
edges. The function W , can be written in terms of polar 
co-ordinates (fig. 2.2) as follows ° 




Fig. 2,2 An edge crack in an infinite plate and 
the co-ordinate system. 
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W = ® + ¥ 

e o 

where w and 'P are even and odd 

o 

function given by 


parts of the stress 


'P = 


Z 


,(-1) 


n-l 


,w 


2 n-l 



O 

cos (n - — ) 8 


2n - 3 
2n + 1 


cos 


(n + 1/2 ) © } 

J 


(2.3) 


+ (~l) n d 2n . r 1 
+ cos (n + 1) 9 ^ j 

and 

n+1/2 

r [sin (n - 3/2 )9 

1 


- sin (n + 1/2) ©] 


If tr,Q) 


o 


n=l, 2 , . .oo 


(- 1 ) 


n-l 


2n- 


f 

1 


- cos (n-l)8 


+ (-l) n a 2n [ - sin (n-l)© + sin (n+1)©] r 

. * . ( 2 . 4 ) 

The constants d n 's and a n 's depend upon. the bouhdaiy 
conditions* more specifically either upon boundary con- 
ditions, at infinity in the case of infinite sector, or 
upon those at some fixed radius when the plate has finite 
dimensions. 


2.5 Formulations of Stress and Displacement Equations! 

Since the "William's stress function is basically 
an Airy's stress function, the three basic stresses can 
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be derived by successive partial differentiation as 
follows : 

d 2 y 3 -r 

; 0. 




3 ^ r 


x 


Dy 2 


y 


a .- 2 


' X 


xy 


3 xB y 

... (2.5) 

In matrix notation these equations can be written in 
terms of polar co-ordinates as follows. 


ri 


"sin 2 © 


sin2©. 

COG © 
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-sin2© 

r 2 

hr' 


, cos 2 © 


~sin2@ 

r 

sm © 
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- 

sin29 
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-sin© 

cos© 

-COs29 

r 

si n© co s© 
r 

cos2Q 
r 2 ' 

where 

M T = 

r c> 2 y 

: L 3 r 2 

d 2 y , 

3 if 'dy 
dr 3© 

32 

3© 

1 ] 

2 J 


r 

?* 

T 


The radial and circumferential displacements, 
V r and U e respectively, can be formulated in terms of 
the biharmonic stress function according to Caber and 
Pilon [s] as follows: 


2 ^ V"§? + Cl - r) r ^ 

*><, n - I SJL 4. (1 v W 2 

V* r cT© + (1 “ y) r 3” 


(2.7) 


where 


ju = ^hear modulus 

r = v/l4v where v is the Poisson's ratio 


and is a harmonic function and is related to the 
stress function W as follows: 


x[0] 


P 

#: 


Making use of the above equations (2*7) and (2.8) the 
displacement field equations, U and U Q can be found out 
easily. To facilitate applications of boundary condi- 
tions U and U_ can be transformed to the cartesian co- 

JO H 

ordinates using the following relationship (fig. 2.3) • 


U = U cos © - U a sin 9 
x r © 


(2.9) 

U s U sin 9 + UU cos 9 
y r © 

where '0' is the angle between the cartesian and polar 
co-ordinate system. 


2.6 Relationship Between Stress Intensity/ Kj , Factor# 

and the First Coefficient of William's Stress 

Function c* 5 
/ "1 


The stress intensity factor , may be derived 
in terms of the first coefficient of the William's stress 
function d^ , as follows: 

From the equation (2.3) it is obvious that the 
* other terms" vanishes at the crack tip i.e. when r * 0. 
Hence# taking only the first term, we get 

CT _ tls- = 'll S£f-2 ( l + sin 

y - g x 2 r 2 2 2 

... ( 2 . 10 ) 


The expression for * <5 y* # qiven in literature based 
on Westergaard [6] is as follows: 
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6 . = 


KL 


/2 7 r 


(1 + 


sin © 


sin 


39 


) ( 2 . 11 ) 


comparing equation (2.10) and (2.11) we find 

Kja - v/~ . d x (2.12) 

This relationship gives us a method that if by some 
means the constants d^ are found out for a particular 
problem then can be determined from d^. The reason 
for not considering the odd series of the stress func- 
tion will be given later along with the method of 
solutions . 


2.7 Single Edge Crack Specimens 


The single edge notch specimen (fig. 2.4a) has 
been analysed by Sraweley et al[3] using the stress 
function ¥/ and its normal derivative &f/3n , as the 
boundary conditions and their results compared well with 
the results obtained by other methods. It is to be noted 
that it was possible for Sraweley to use the derived 
boundary conditions ('P and ^ ?/5 n) as the loading and 
geometry of the specimen were fairly simple. The boun- 
dary conditions used by Srawley are as follows (fig. 2. 4a) 
Along A-B 


¥ = 0 } 
Along B-C 


5_f 

2 > x 



¥ = <r„ ( 


x + a\2 


o 


d ¥ 
2>y 


= 0 ; 


(2.13) 
contd. . . 
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where a, and W are the crack length and width of the 
plate and 5"" is the uniform tensile stress applied 
along B-C. 

In our work-in an attempt to use the traction, 
and displacement boundary conditions involving second 
order partial diff erentiation ; the following direct boun- 
dary conditions are used. 


Along A-B 

as 0 ; Z = 0 ; 

w x ' xy ' 

Along B-C 

6 = 6 i £• t= 0 ; (2.14) 

y o' u xy 

Along C-D 

6 ' at 0 s ?" = 0 s 

x * ‘-xy 

As can be observed from the figure there exists a symmetry 
with respect to the x-axis and hence only top half the 
specimen and hence the even part of the William's stress 
function need to be considered. Therefore 

'i' (r, 9) a ? (r, ©) (2.15) 

This is true for all the cases where there is symmetry 
about x-axis and since in our present work we deal only 
problems of this kmd^ eguation (2 .15 ) is true for all 
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the cases that we consider. Also if <35 ^placement condi 
tions are prescribed along the boundary B-C (fig. 2.4b) 
then the boundary conditions become as follows: 

Along A-B 




Along B-C 


U x = 0 * 


Along C-D 


6 x ° 0 ' 


T = 0s 
°xy ' 


U = v ; 
y o 


t = 0 ; 

xy ' 


(2.16) 


A square specimen of this kind is analysed by Kobayashi 
[ 7 ] using finite element techniques. 


2.8 Centre Crack Problem: 

The centre crack problem is treated as two edg „ 
crack problems which are separated by the mid plane L~T/ 
(fig. 2.5a) as shown. We consider only one quarter of 
the specimen because of symmetry, with the following 
boundary conditions (fig. 2.5b) 

Along A~B 

u x = 0 > £xy = 0 » 

Along B-C 

^y- 6 "o> txy = °> (2 - 17) 

Along C-D 

^=0 ; t xy . 0 



Ul'l 



Rg. 2,5 Co) Centre crock problem (b) One 'quarter 
©I the specimen tor analysis purpose. 
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This specimen is analysed by M« Isida and others [ 2 ] 
using a complex stress function prescribed for centre 
cracks. In our work we use the traction and displace- 
ment boundary conditions directly instead of ¥ and ^/£>n 
as used by Isida. We use the same William's stress 
function which is prescribed for edge cracks, for this 
centre crack problem incorporating the displacement boun- 
dary condition. Thus the usage of William's stress func- 
tion is extended to the centre crack problems as well in 
our work. 

2.9 Problem of Butt-Weld with Improper Penetration: 

Weld defects constitute a group of stress raisers 
which have more serious effects than reinforcement or 
undercuts, of the weld profile. One category of such 
defects covers those caused by insufficient penetration 
in butt-welds. The American welding society [8] reco- 
mmends that the butt-welds with partial joint prepara- 
tion be used for static loading only. According to their 

tests full penetration butt-welds on plates (square 

"tuh. 

ground welds) are only possible on plates 1/8 of an 
inch. With a root opening of one half the plate thick- 
ness, full penetration is obtainable upto 1/4 inch thick, 
but the root of the first weld must be chipped out to 
sound metal before depositing the second weld. For 

these recommendations can not 


economical reasons. 
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always be adhered to and butt welds with improper pene- 
tration are common in non-critical welded connections 
subjected to cyclic loads. The economical need is felt 
for standards that specify permissible defect sizes in 
welded metal structures. Lack of penetration is assumed 
to be similar to a crack at the interface of two plates 
joined by welding (fig. 2.6a). Using the stress inten- 
sity factor for this geometry and loading a very good 
approximation can be made for the stress field and then 
we can determine whether a joint will perform to our 
expectations • 

The boundary conditions for this geometry (fig. 
2.6b) is as follows! 



Along DEF 

<T n - 0 » - 0 » 

where 5~ n and £ g are normal and shear stresses along 
the curved boundary which is circular with radius R 
DEF, and are related to (T x * CT y ^ Z xy stresses as 



&« 



Fig. (a) Butt-weld improper penetration problem. . 

(b) . One quarter of the butt -weld problem for 

analysle purpoee. 

(c) Co-ordinate repreeentatlon of etrees transfer-, 
matione for the curved boundary *OEF' 



follows (fig. 2.6c) 
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CX n = o~. cos 2 a + 2<^ y sina cosa + sin 2 a 

r s = -° y ) sina cosa - ~C ^ (cos 2 a - sin 2 a) 

... 2 . 19 

and a is as shown in (fig. 2 f 6c). 

We discuss the method of solution and the 
results for all the problems formulated^in the next 
chapter* 



CHAPTER - 3 


RESULTS AND DISCUSSION 

3.1 Introductions 

In this chapter we describe the method of solu- 
tion and the results in detail as applied to single 
edge-crack problems and centre crack problems. The 
results obtained agreed well with the existing results 
only upto H/W ratio of 0.6, thus necessitating the need 
to improve the solution procedure. An improvement in 
the solution procedure is suggested in the present work, 
and is applied to all the problems discussed in sections 
(2.7, (2.8) and (2.9). Reasonably good results were 
obtained even upto H/W ratio of 1*5. The results are 
presented in a normalized form. 

3.2 Solution Procedure and Discussion of Results: 

If we consider the single edge notch problems 
and centre crack problems as discussed in sections (2.7), 
(2.8) and (2.9), we observe that, at each point over 
the boundary of the specimens we have two boundary _ 
conditions to be satisfied as given by equations 2.14, 
2.16, 2.17 and 2.18. 
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Since we have to take only the even series of the stress 
function, , as explained earlier, because of symmetry, 
we get exactly '2n' constants for each 'n', in the res- 
pective expressions of the boundary conditions. Thus 
to solve the problem we truncate the expressions with 
* 2n' arbitrary constants (d n 's) and satisfy the boun- 
dary conditions at 'n* number of points selected over 
the boundary at equal distances. We get exactly # 2n' 
equations in / 2n / constants. Nov; if we solve this sys- 
tem of linear equations we can find stress intensity 
factor, K_ , using the relation 2.12. 

3.2.1 Selection of Boundary Points: 

A subroutine is developed which would give the 
co-ordinates of the required 'n' number of points over 
the boundary say A.B.C.D, in the case of a single-edge 
notch problem (fig, 2.4a) distributed at equal distances 
along the periphery. The input parameters to this 
routine are, the crack length, a, width of the plate 
'W', the height of the plate from the crack plane H, 
and of course, the required number-of points, *n ' . Thus 
the geometry of the plate is taken care of the co-ordi- 
nates of the boundary points that subsequently enter 
main solution procedure* It can be stated tyere that 
the input preparation to the < 0110 cation p^6cedure is 
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very simple when compared to the finite element method. 
3.2 1 , 2' Solution of Equations: f 

The boundary conditions are satisfied exactly 
at the ' n ' number of points, i.e. the error is made zero 
at these points in Collocation procedure. Gaussian 
elimination is employed to solve the system of linear 
equations. We are interested in finding mainly, the 
first coefficient of the William's stress function d^, 
which is to be used to compute the stress intensity 

I 

factor, Kj» If we bring the constant d^ to the last 
place of the truncated expressions Gaussian elimina- 
tion eliminates the possibility of additional computa- 
tions and thereby subsequent error in finding d^. 

As we are approximating the entire boundary 
points by a finite number of boundary points, n, it is 
expected that the exact results would be approached 
when , n t tends to infinity. Also we approximate the 
stress field by using , ’k t number of constants such 
that k ss 2n. Reasonably good results are obtained 
within allowable error if the number of points * n* 
increased until the coefficient d 1 reaches a stable 
value. This procedure is repeated for differqnt values 
of for a particular plate dimensions i.e. fixing 
a/W and H/W. The solutions for single edge crack 
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problem is presented in a tabular form in which the 
value of d^ is shown against for each f n* number of 
points taken over the boundary (table 3. a). 

It is observed that this method of solution 
i.e. taking exactly the required number of points over 
the boundary to find the unknown coefficients d's does 
not give good convergence of when the ratio H/W > 

0.6; (table 3. a) therefore the need to improve the 
technique is felt and the following method is suggested. 

We truncate the expressions of the boundary 
conditions to *k' number of constants and select * n' 
number of points over the boundary such that 2n > K. 

That is we take more number of points than required to 
find the constants and treat the system of linear equa- 
tions as over determinant. This kind of treatment of 
the system gives good convergence of the coefficient 
d 1 as both the number of points ' n * and constants 'k' 
are increased. The values of for different a/W and 
H/W ratios against the number of boundary points 9 n ' 
and arbitrary constants d* are presented in table ( 3.b) • 
We employed subroutine 'bLSQAR' available in our 'IMSL' 
routines. This method of solution gives good results 
even upto H/W ratio of 1*5, above which the stress inten- 
sity factor is assumed to be insensitive for all prac- 
tical purposes as there is no appreciable change observed. 
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TABLa-?a 

Results obtained by solving (2n x 2n) system for 
single- edge-crack problem with traction prescribed 
boundary conditions 


Trial a/W 

No. 

H/W 

Nq , of 
points 
over the 
boundary 

n 

No. Of 
cons- 
tants 

'k' 

*Mrst co- 
efficient 

d l 

1 

—"5 

3 

j 


6 

i 

" n o.S‘ 

6.4 

1 18 



16^3* 

2 

0,5 

0.4 

21 

42 

*•* 162 * 62 

3 

0.5 

0.4 

24 

48 

- 163.87 

4 

0.5 

0.4 

27 

54 

- 165.31 

5 

0.5 

0.4 

30 

60 

- 164.95 

i 

0.5 

0.5 

18 

36 

- 167.23 

2 

0.5 

0.5 

21 

42 

- 134.1091 

3 

0.5 

0.5 

24 

48 

- 150.56- 

4 

0.5 

0.5 

27 

54 

- 150.9178 

5 

0.5 

0,5 

30 

60 

- 150.3091 

1 

0.5 

0.6 

la 

36 

- 140.96 

2 

0.5 

0.6 

21 

42 

- 125.21 

3 

0.5 

0.6 

24 

48 

- 146.075 

4 

0.5 

0.6 

27 

54 

- 145.29- 

5 

0.5 

0.6 

30 

60 

- 147.847 

1 

0.4 

0.6 

18 

36 

- 94.213 

2 

0,4 

0.6 

21 

42 

- 76.222 

3 

0,4 

0.6 

24 

48 

- 99.15 

4 

0,4 

0.6 

27 

54 

- 99 .02 

5 

" 0.4 

0.6 

30 

60 

- 101.43 
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TABLx,-3b 


Results obtained by solving (2n x k) system for 
single edge crack problem with traction prescribed 
boundary conditions 


Trial 

a/W 

H/W 

No. of 

No. Of 

First coeffi- 

i A * i tl 


. / . S5- 

points 

cons- 

cient 



"V > 

A 

over the 

tants 




X 

boundary 


dL 



\ 

n 

k 

1 

’ 1 

2 

3 

4 

5 

6 




75 

50 

- 63.216 




73 

50 

— 62.32 3 




80 

50 

— 63* 140 


0*3 

1 .0 

83 

60 

- 63.165 




80 

60 

- 62.99 




83 

60 

- 62 . 300 




30 

50 

- 93.96 




83 

50 

- 93.93 


0.1 

1*0 

80 

56 

- 94.13 




83 

56 

- 94.14 




80 

60 

- 141.26 


0.5 

1.0 

83 

60 

- 141.06 




80 

60 

- 219.09 


0.6 

1,0 

83 

60 

- 220,99 




80 

56 

- 258,91 




83 

56 

- 376.82 




80 

46 

- 314.75 


0.7 

1.0 

83 

46 

- 315.82 




80 

44 

— 295.02 




33 

44 

- 292,28 




7. 

50 

— 64.75 



- 

n 

50 

- 64.88 




% 

60 

- 64.89 


0,3 

0.5 

78 

60 

— 64.948 




80 

60 

— 64.91 




83 

60 

- - 64.80 

- 

■ 


75 

50 

- 95.340 


0.4 

0,5 

78 

50 

- 95.322 






Oontd. 



Contd 


.Table- 3b 
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2 

3 

4 

0.5 

0.75 

75 

7 s 

0.6 

0.75 

75 

78 

0.7 

0.75 

75 

78 

80 

83 

80 

83 

80 

83 

0.3 

0.5 

75 

78 

75 

78 

75 

78 

0*4 

0.5 

75 

78 

80 

83 

0.5 

0.5 

80 

83 

80 

0.6 

0.5 

75 

78 

0.7 

0.5 

75 

78 

0.2 

0.5 

75 

78 


5 6 


50 


142.04 

50 

*•* 

142.04 

50 

«w 

221.25 

50 

mm 

221.25 

50 

mm 

374.51 

50 

mm 

368.57 

50 

mm 

375.6 

50 

mm 

368.89 

60 


377.05 

60 

mm 

370.07 

46 


369.70 

46 

- 

369.163 

40 

- 

70.96 

40 

mm 

71.19 

50 

mm 

71.52 

SO 

mm 

71.43 

60 

mm 

71.63 

60 

- 

71.51 

60 

mm 

104.02 

60 

- 

103.93 

60 

m 

103.90 

60 

mm 

103.95 

60 

mm 

150.40 

60 

mm 

150.42 

64 

mm 

150.41 

60 

mm 

227.19 

60 

mm 

227.11 

60 

mm 

377. sa 

60 

mm 

377.99 

60 

mm 

46,68 

60 - 

- 

46.52 




' J S 






^ dteti » w js® 
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stress intensity factor is presented as a 
dimension! ess ratio against a normalized crack 

length a/W, fig. (3.1), The normalizing parameter Kg 
has the dimensions of a stress intensity factor and is 
so chosen as to emphasize some important aspect of the 
solution, e.g, to emphasize the effect of additional 
boundaries, would be equal to K 3 in the absence of 
such boundaries [ 9 ]. This fgrm of presentation shows 
clearly the effects of important features such as geo- 
metrical boundaries and stress distributions. 

The improved method of solution is applied to 
the problem discussed in sections 2,7 to 2.9 and their 
results are plotted in the normalized form as shown in 
figs, (3,2), (3,3) and (3,4). The convergence of the 
coefficient d^ against number of boundary points # n # 
and the number of constants # k' for different a/W and 
H/W ratios are also presented in tables (3,c) (3,d) and 
(3.e). A square specimen of single edge crack with 
di splacsement prescribed boundary conditions (section 
2,7) has been analysed by Kobayashi [#?] et al # using 
finite element method# In our present we studied the 
specimen with three different H/w ratios of 0,5, 0,75 
and 1*0, For H/w ratio of 0.5 our results compared 
well with those obtained by finite element method for 
a/W >0*2. The maximum error limit when compared with 
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TABLE- 3c 

Solution for single edge notch specimens with 
displacement prescribed boundary conditions 

H/W « O.S ; V Q s 100.0 


a/w 

NO. o€ 

boundary 

points 

'n' 

NO • Of 
constants 

'k' 


~T 

& 

V i 

4 


30 

30 

- 30.13 

0.2 

33 

30 

- 29.82 


30 

35 

- 29.88 

0.3 

33 

35 

- 30.10 

* 

30 

30 

- 30.618 

<** 


33 

30 

- 30.519 

0.4 

30 

40 

- 31*566 


33 

40 

- 30.350 


30 

35 

- 31.520 

0,„5 

33 

35 

- 31.874 


30 

36 

- 32.21 

0.6 

33 

36 

— 32 . 34 

- 

30 

35 

- 32.93 

0.7 

33 

35 

- 32.89 


Contd. , . .Table- 3c 


H/W = 0.75 


1 

2 

3 

4 


30 

39 

- 26.816 

0.2 

33 

39 

** 2 6 *88 6 


30 

39 

- 28,058 

0.3 

33 

39 

- 28.112' 


30 

30 

- 28.917 

0.4 

33 

39 

- 29.956 


30 

39 

- 29.558 

0.5 

33 

39 

- 29.570 


30 

39 

- 29,743 

0.6 

33 

39 

- 29,668 

- 

30 

39 

- 29.726 

0.7 

33 

39 

- 29.534 


30 

39 

- 31.133 


33 

39 

- 29.891 

0.9 

30 

35 

- 30. 96 a 


33 

35 

- 30.425 


Contd.... 


39 


Contc3« . 

. .Table- 3c 

H/W « 

1.0 


1 

2 


3 

4 


30 


35 

- 23.035 

0.2 

33 


35 

- 2 3.17 3 


30 


35 

- 24.985 


33 


35 

- 25.144 

0.3 

40 


49 

- 25.177 


43 


49 

- 25.326 


30 


39 

- 26.055 

0.4 

33 


39 

- 26.212' 


30 


39 

- 26.700 

0.5 

33 


39 

- 26.834 


30 


39 

- 26,951 

0.6 

33 


39 

- 27.050 

- 

30 


39 

- 27.114 

0.7 

33 


39 

- 27.197 

0.8 

30 


39 

- 27.987 


33 


39 

- 27.3 36 

0.9 

30 


39 

- 28.959 


33 


39 

- 28.534 


40 


TABLE- 3d 


Solution for centre crack problem with traction 
prescribed boundary conditions 


« 100,0 ksc. ; HA* * 0,5 

V 


a/W 

( 1 ) 

No. of boundary 
points, n 

(7) 

No. Of 
constants 

k (.3) — 

d i 

ft) 


30 

50 

- 23*822 

u e x 

33 

50 

- 2 3.692 


30 

48 

- 37*151 

0.2 

33 

48 

- 37.128 


30 

48 

- 5 3.091T 

0,3 

33 

48 

- 53.091 


30 

50 

- 72.872 

0.4 

33 

50 

- 72.868 


30 

50 

- 98.375 

0.5 

33 

50 

- 93,361 


30 

50 

-132.826 

0.6 

33 

50 

-132.787 


30 

50 

-179.699 

0.7 

33 

50 

- 179.799 


30 

50 

-217.2755 


33 

50 

-223.683 

0.8 

30 

20 

-243.53 


33 

20 

-2 30 ,50 


Gontd* » * • 


Contd. • •• .Tab! e~ 3d 


H/W » 0.7 

' <^0 

=. 100.0 

kse 

1 

2 

3 

4 


25 

40 

- 23.921 

0.1 

28 

40 

- 23.892 


35 

60 

- 24.041 


38 

60 

- 24.048 


35 

60 

- 34.875 

0.2 

33 

60 

- 34.374 


35 

60 

- 47.55 4 

0.3 


60 

60 

- 47.555 


35 

60 

- 62.601 

0.4 


38 

60 

- 62.602 


35 

60 

- 30.941 

0.5 

33 

60 

- 30.940 

o.g 

35 

60 

-103,164 

33 

60 

-10 3.161 


35 

60 

-129.483 


38 

60 

-129,896 

0,7 

40 

64 

-126.761 


43 

64 

-129,202 


Contd, .. . 


Contd Table 3d 


H/W =1*0 , cr 0 = 100.0 *SC 


1 

2 

3 

4 

o.l 

30 

50 

— 2 3*8 32 


33 

50 

- 24.059 


30 

50 

- 33.487 


33 

50 

- 33.531 

0 .2 

35 

50 

— 33.456 


38 

50 

- 33.426 


35 

60 

- 33.042 


33 

60 

- 33.570 


30 

48 

- 43.511 

0.3 

33 

48 

- 43.516 

0.4 

30 

48 

- 54.386 

33 

43 

- 54 . 332 


30 

43 

- 66.687 

0.5 

33 

48 

• 66.633 


30 

48 

- 81.074 

0.6 

33 

48 

- 81.098 


25 

40 

- 98.923 


28 

40 

- 99.102 

0.7 

30 

48 

-101.722 


33 

43 

- 99.923 


Contd 
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Contdo o . . .Tab! e~3<3 


H/W=l .5 $ ff Q « 100.0 ksc 


1 

2 

3 

4 


35 

40 

- 26.975 


33 

40 

- 25.090 


45 

30 

- 26,714 

0.1 

48 

30 

- 26.690 


40 

30 

- 26.688 


43 

30 

- 26.730 


40 

40 

- 33.061 

0 0 2 

43 

40 

- 30 . 561 

40 

30 

- 34.209 


43 

30 

- 34.210 


40 

30 

- 41.939 

0 „ 3 

43 

30 

- 41.929 


40 

30 

- 50.404 

^ .4 

43 

30 

- 50.380 


40 

30 

- 60.071 


43 

30 

- 60.0 35 

to 

c 

o 

40 

40 

- 60.060 


43 

40 

- 60.075 


40 

40 

- 70.089 

0.6 

43 

40 

_ 70.307 


40 

40 

- 82.786 


43 

40 

- 92.564 

0.7 

<10 

30 

- 84.514 


43 

30 

- 84.587 


40 

34 

- 36.224 


43 

34 

- 86.173 


40 

36 

- 86.907 
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TABLE 1 * 3e 

Solution, for Butt weld improper penetration problem 
W e 0.5 inch 
R a 1.04 inch 
C a 0 . 3 inch 


a/w 

% 

"V 

No .of boundary 
pblnts 'n' 

No. Of 
constants 

'k' 

first 

coefficient 

d l 


25 

40 

*24,203 

0.2 

23 

40 

*24.180 


25 

40 

*30.993 

0.3 

23 

40 

*30.991 


25 

40 

*33.410 

0*4 

23 

40 

* 19*410 


25 

40 

*46.571 

0.5 

23 

40 

*46.572 


25 

40 

-55.709 
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finite element method for this range of solution is 
4*5#. The solution of this problem is plotted in fig* 
(2.2) with finite element solution. 

The centre crack problem under tensile loading 
as discussed in (section 2.8) is analysed and the results 
(fig* 3.3) compared well with those obtained by M.Isida 
[ 2 ]. M, Isida used a complex stress function prescribed 
for centre cracks using the stress function and its 
normal derivative as the boundary conditions. The good 
agreement of our results with those of Isida' s gives 
us confidence to proceed further with the problem of 
butt-weld partial penetration. The result of this 
problem for a sample input of H/W ratio equal to 1.0 is 
plotted in fig. (3.4) in the normalized form. The solu- 
tion of stress intensity factor can be used to find the 
stress concentration factor, by assuming a fictitious 
crack tip radius as follows: 

K x o i Kjj. 6' . /Tf (3.1) 

where is the elastic stress concentration factor 
and r is the far field' stress from the crack tip. 

The stress intensity factor K^. can also be used to study 
the fatigue growth of the cracks and permissible defect 
sizes can $lso he predicted. 
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3*3 Conclusions: 

From our present work on Collocation procedure 
to find stress Intensity factors for both single edge 
crack and centre crack problems the following ©onclu— 
sions can be made. 

1) Collocation method with boundary conditions 
involving second order partial differentiation of the 
stress function give good results for the problems dis«- 
cussed in on* work,if the number of boundary points 
taken is more than the number of arbitrary contents of 
the William's stress function. 

2) Collocation procedure offers the distinct advan- 
tage of easy input preparation. For different a/w and 
H/W ratios, solutions are obtained without tedious input 
preparations. 

3) The usage of William's stress function is exten- 
ded to centre crack problems as a step to generalize 

the procedure. 

4) The convergence of the coefficient# d^ is good, 
as the number of points over the boundary and the number 
of constants of the stress function are increased even 
if a part of geometrical boundary is curved (e.g. weld 
specimen) • 
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5) This method can be used to optimize the speci- 
men size for fracture toughness testing purposes. For 
example in the case of single edge notched specimen 

as there is no appreciable change in stress intensity 
factor for ratio of H/W above 0.75 for a/W ratios 0.2 
(fig. 3.1) the specimen size can be limited with H/W 
ratio as 0.75. 

6) The easy method of finding stress intensity 
factor for weld specimens enables one to carry out 
studies on fatigue growth of cracks and to prescribe 
allowable defect sizes of weld geometries. 
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